APS/123-QED 



Asymmetric Two-dimensional Magnetic Lattices for Ultracold Atoms Trapping and 

Confinement 

A. AbdelrahmanQ M. Vasiliev, and K. Alameh 

Electron Science Research Institute, Edith Cowan University, 
270 Joondalup Drive, Joondalup WA 6027 Australia 

P. Hannaford3 

Centre for Atom Optics and Ultrafast Spectroscopy, 
and ARC Centre of Excellence for Quantum Atom Optics, 
Swinburne University of Technology, Melbourne, Australia 3122 
(Dated: October 27, 2009) 

A new method to implement an asymmetrical two-dimensional magnetic lattice is proposed. The 
asymmetrical two-dimensional magnetic lattice can be created by periodically distributing magnetic 
minima across the surface of magnetic thin film where the periodicity can be achieved by milling 
n x n square holes on the surface of the film. The quantum device is proposed for trapping and 
confining ultracold atoms and quantum degenerate gases prepared in the low magnetic field seeking- 
state at low temperature, such as the Bose-Einstein Condensate (BEC) and ultracold fermions. We 
present detailed analysis of the analytical expressions and the numerical simulation procedure used 
to calculate the external magnetic field. We also, describe the magnetic band gap structure exhibited 
by the asymmetric effect of the magnetic minima and show some of the possible application. We 
analyze the effect of changing the characteristic parameters of the magnetic lattice, such as the 
separating periodicity length and the hole size along with the applications of the external magnetic 
bias fields to maintain and allocate a suitable non-zero magnetic local minima at effective z-distance 
above the thin film surface. Suitable values are shown which keep the trapped ultracold atoms away 
from the thermal Majorana spin-flip and the surface Casimir- Polder effect. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

The emerging field of quantum computing and infor- 
mation processing using ultracold atoms prepared in spe- 
cific quantum states and trapped in specifically patterned 
magnetic minima, has attracted a great interest over the 
last decade where magnetic micro-traps have been re- 
cently recognized to have more applicable accessibility 
to several degrees of freedom of their individual hosted 
quantum systems [? ]. Several micrometer-scale struc- 
tures generating magnetic periodic potential have been 
proposed for such purpose as an alternative to optical 
lattices [? ]. It is possible to produce a lattice con- 
figuration by periodically introducing specific patterns 
on the surface of a magnetic material, where in such 
configuration they are known as magnetic lattices [? ? 
? ? ? ], in which they can be realized by manufac- 
turing one-dimensional and two-dimensional periodically 
distributed magnetic minima. Such specifically engi- 
neered quantum devices are also classified as Atom Chips 
[???]. Magnetic lattices can also be created using 
current carrying wire, however a few technical difficul- 
ties have been encountered in relation to the integration 
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of current-carrying wires onto the magnetic micro-traps. 
For example, current-produced near-field thermal noise 
reduces the number of trapped cold atoms. Also, the 
current density fluctuations lead to the fragmentation in 
the BEC cloud [? ? ]. These technical difficulties may 
be avoided by using optimized magnetic materials pos- 
sessing high remanent magnetization and enabling highly 
stable trap frequencies through a tight magnetic confine- 
ment and reproducible periodic potentials [? ]. It may 
also offer a controllable phase transitions, possible ma- 
nipulation and less demanding access to the quantum 
ensemble for measurement [? ]. Also, the ability of mag- 
netic lattices to attain pronounced confinement and a 
proper control of single particles makes them excellent 
candidates to host large-scale qubit systems [? ? ? ? 
]. In this article, we propose, a simple method to imple- 
ment an asymmetric two-dimensional magnetic lattice, 
using the current micro- fabrication technologies, where 
we analyze the performance of the magnetic lattice for 
ultracold atoms confinement and trapping, with the aid 
of applying in situ external magnetic bias fields that are 
integrated on the atom chip. In section (|TT]) we derive 
the analytical expressions which describe the magnetic 
minima distribution across the surface of the magnetic 
thin film. In section (|III[) we show the simulation proce- 
dure used to calculate the surface effective magnetic field 
and in section (|IV[) we describe the relevant characteris- 
tic parameters that will be used to optimize the magnetic 
minima distribution and hence preparing suitable quan- 
tum states such as two-level systems. We conclude in 
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section (jVj by showing the possible application of this 
type of magnetic lattice. 



II. MAGNETIC LATTICES STRUCTURE 

The proposed structure to generate the magnetic lat- 
tice based quantum device is realized by milling an m x m 
array of blocks, each block being an array of n x n square 
holes of size ah , separated by a s , where n represents the 
number of holes in a block as detailed in Figures lll.lf a- 
b). The gaps between the blocks containing no holes 
are assumed to be greater than, or equal in width to 
a s . The holes are formed, within a magneto-optic thin 
film of thickness Tbtm sputtered onto a proper substrate. 
The depths of all holes are equal and extend through the 
thin film thickness down to the substrate surface level. 
The gaps separating different blocks are an important 
design feature which enables the n x n magnetic lattices 
to be surrounded by an unperturbed film area, which 
introduces an extra degree of confinement through the 
creation of magnetic field "walls" encircling the n X n 
matrices, and isolating them from one another as shown 
in Figure Hi. lf hV The magnetic quantum device is in its 
remanently-magnetized state, where its magnetization, 
M z , is directed in a direction perpendicular to the sur- 
face plane. Effective parameters, i.e., a/,, a s and rum of 
the magnetic lattice for different values of n are consid- 
ered. 

The structure generates two dimensional periodically 
distributed magnetic field minima in the vicinity of the 
surface of the perforated film where the distribution cre- 
ates the magnetic lattice that used to trap the cold atoms. 
The presence of holes results in a magnetic field distribu- 
tion whose non-zero local minima are located at effective 
z-distances from the holes opening centers above the film 



plane, denoted here by d m i n . These minima are localized 
in small volumes representing the potential wells that 
contain certain number of quantized energy levels for the 
cold atoms to occupy. In our design, we assumed that 
the size of the holes ah and the holes separation a s are 
both equal, ah — a s = a, to simplify the mathematical 
derivation where we adopted an analysis approach simi- 
lar to that reported in [? ? ]. The spatial magnetic field 
components B x , B y and B z can be written analytically as 
a combination of a field decaying away from the surface of 
the trap in the z-direction and a periodically distributed 
magnetic field in the x — y plane produced only by the 
magnetic induction, B a = [i M z /tt, at the surface of the 
permanently magnetized thin film. We define the surface 
reference magnetic field as B re f — B Q (1 — e~^ r ), where 
f3 = n/a, t = Tbtm is the film thickness, and the plane of 
symmetry is assumed at z = 0. Although our model 
takes into account the magnetic local minima created 
only by B Q (similar to a flat magnetic mirror), we also 
include in the analysis the external magnetic bias field 
components along x, y and z directions, namely, B x -bias, 
By-bias and B z - ias which can also be localized by means 
of micro-fabrications. Figure HLir d-el- TH.lf f-g) show the 
simulated map of the magnetic field strength distribu- 
tion across the x — y planes located at different distances 
above the magnetized film surface. Both the density plot 
and contour plot representation types are shown for both 
simulation results of a 10 x 10 magnetic lattice at the ini- 
tial magnetic state, formed by B re f only, and with the 
external application of B x - ias and B y - ias, respectively. 
The analytical expressions that describe the non-zero lo- 
cal minima, periodically positioned in the x — y plan, 
take into account the strength of the effective field and a 
[? ]. These expression are derived and simplified to the 
following set of equations 
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The higher order terms in these equations, which describe 
the far field, can be neglected for distances that are large 
compared to the effective distance. This is because only 
the cold atoms that are in the so called low magnetic field 
seeking-state (atom's magnetic moment to be oriented 



antiparallel to the localized magnetic field in the trap) are 
attracted to the non-zero local minima of the magnetic 
traps located at the effective z-distance from the film 
surface, which is larger than a/2ir . Thus, Equations (JTJ) 
to are reduced to the fallowings 
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FIG. II. 1: (a) A 10 x 10 magnetic lattice surrounded by an unperturbed area, (b) The lattice parameters are specified by the 
hole size oth x an, the periodic spacing a s between the holes and the magnetic layer thickness Tbtm- (c) Magnetic density plot 
of the simulated magnetic lattice sites at z — x plane along the center of the lattice. The traps (dark color) are located at 
the effective z-distance, d m in, from the holes opening centers of high magnetic field (bright color) above the film plane, (d-e) 
Magnetic field density plot across the x — y plane at the d m in with no external magnetic bias fields and with application of 
the external B x -t, ias — 10 G and B y -t,i as — 10 G magnetic bias fields, respectively, (f-g) Contour plot of distributed lattice 
sites across the x — y plane without and with applications of B x -tias — 10 G and B y -t,i as — 10 G respectively, (h) 3D plot of 
the magnetic field (on the z-axis) of the distributed sites across the x — y plane (a;-axis and y-axis) simulated at the dmin and 
displayed from the center sites to the edge sites where B x -i, ias = 10 G and B y -uas — 10 G were applied. Simulation input: 
a s = a h = 1 [im M z = 3 kG, r btm = 2 \im and T p - wa ii = 1 fim. 



B x — 



B„ — 



x — bias 



B (l - e'^j e-^-^smiPx) + B 
B (\ - e^J e-^-^sin(/3y) + B v . b 



B z = flJl-e-Me-^] 



The magnitude B of the magnetic field above the film 
surface can be written, using Equations (HJ) to ©, as 
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The simulation result of numerical calculations for a fi- 
nite magnetic lattice are shown in Figure III. 21 Figure 
III.2f a) shows the location of the magnetic field non-zero 
local minima along the z-axis, B^ in , at the effective dis- 
tance ,d m i n , from the hole opening center and confined 
along z-axis by a magnetic barrier of magnitude AB Z . 
The results shown in Figure III. 2T b) demonstrate the ex- 
istence of non-zero local minima of the magnetic field 
at the effective z-distance along the x-axis, confirming 
the suitability of the structure proposed for trapping the 
atoms being in the low magnetic field seeking state. The 
intra-block magnetic lattices consist of tunneling barri- 
ers separating the trapping sites where each magnetic 
quantum well characterized by the non-zero local min- 
ima and differs from each other by amount of tilted po- 
tential AV as shown in Figure llL2T d). The unperturbed 
film area, around the n X n hole-matrix block, produces 
magnetic walls that localize and magnetically decouple 
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the individual blocks from each other, as illustrated in 
Figure HL2T c) . This allows individual quantum-state ad- 
dressing for each nxn magnetic lattice separately, where 
different quantum features can be investigated, such as 
quantum multi-particle (multipartite) long-range entan- 
glement and strong correlations [? ], the large scale lat- 
tice spacing in our magnetic trap is suitable for single 
atoms detection [? ] . It is also possible to vary the indi- 
vidual tunneling barrier heights by applying an external 
magnetic field in the negative z-direction. The numerical 
calculations in Figure llL~3T c-d) show the effect of applying 
small external B z -bias magnetic fields on the tunneling 
barrier across the lattice sites. The z-bias field is possible 
to be applied uniformly from outside the vacuum system 
of the BEC experiment. It is also possible in case of large 
a s to control the tunneling barriers in situ through micro- 
coils surrounding each hole. Table (|T| summarize the 
important parameters used in simulating the proposed 
magnetic lattice-based quantum device, where Tbtm-, T waii 
represent the bottom layer and the wall thickness of the 
two dimensional magnetic lattice, respectively. Table ([IT)) 
shows the output of the two-dimensional magnetic lattice 
characteristic parameters describing the suitability of us- 
ing the quantum device to provide a longer lifetime for 
the trapped ultracold atoms to process quantum infor- 
mation. 
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FIG. II. 2: (a) The magnetic field minima B min along the z-axis at the effective z-distance, d m i„, away from the surface of the 
thin film, (b) The field minima calculated along z-axis at the center and the edge lattice sites, (c) Magnetic field distribution 
simulated along the x-axis for a 20 x 20 two-dimensional magnetic lattice, (d) Two magnetic quantum wells at the lattice edge 
along the y-axis separated by the magnetic tunneling barrier AB and have different values of magnetic non-zero local minima 
through tilted (magnetic) potential AV. Simulation input of the magnetic lattice characteristic parameters: M z = 3 kG and 
Tbtm — 2 pjm with no external bias fields applied. 



TABLE I: Simulation parameters used to produce the non-zero local nil 111 111 £L Bmin with their effective distances d m in- The 
magnetic walls are produced by the condition r wa ii Tbtm- 

Simulation parameters m x m n x n M z Tbtm T wa u a s 

Input values lxl 11 x 11 2 kG 2 fim 1 fim 1 fim 1 /im 



III. SIMULATION OF THE SURFACE 
MAGNETIC FIELD OF THE MAGNETIC 
LATTICES 

To simulate the magnetic lattices, field-producing vol- 
umes are represented by symmetrical segments, where 
each segment is subdivided into a suitable number of 
smaller objects for accurate computation. Uniform mag- 
netization M z is assumed to exist within each segment 
of the total volume of the n x n magnetic lattices, which 
produces the magnetic field distribution in the space sur- 
rounding the structure. Mutual interactions between the 
segmented volumes are accounted for, and all relevant 
data are stored in a relatively large matrix. The simula- 
tion is carried out through the relaxation procedure by 
first loading the material properties, the magnetization 
M z , the layer thickness r = Tbtm, an d the structure ge- 
ometry which is done by specifying the values of ah, a s , 
and n. The external surface magnetic field is then calcu- 
lated at different points in the space around the surface 
of the magnetized thin film. The computed fields at the 
non-zero local minima then provide accurate estimates 
for the depths of the potential wells and their locations 
at the effective z-distance. The numerical results are cal- 
culated using the package RADIA [? ]. 

We simulate the important effects of changing the 
structure dimensions on the magnetic lattice characteris- 
tic parameters, in particular the location of the non-zero 
local minima, d rn in, from the surface of the thin film, 
their sizes and the curvature of the magnetic field across 
each individual magnetic potential well that forms a lat- 
tice site. The gradient of the magnetic field can be con- 
trolled via change in ah , also shallow or steeper magnetic 
potential wells can be realized when choosing the suit- 
able dimensions of a s and ah- Figure (|H.3|) summarize 
the simulation results of changing the hole size ah, the 
length of the periodicity a s and applications of the exter- 
nal magnetic bias field along the negative direction of the 
z-axis, B- z -bias, and their effects on the relevant charac- 
teristic parameters including the height of the tunneling 
magnetic barrier AB and the amount of the potential tilt 
AV between two adjacent sites of the magnetic lattice. 



IV. CHARACTERISTIC PARAMETERS OF 
THE ASYMMETRICAL TWO-DIMENSIONAL 
MAGNETIC LATTICE 

The initialization setup of the magnetic lattice, using 
only the magnetic induction at the surface of the per- 
manently magnetized thin film in the z direction, leads 
to a symmetrical distribution of the potential wells in 
the x — y plane. We simulate this effect by applying 
the constrains B x ^ bias = B y - blas = B z - bias = and 
B ~ /i M z /tt, where the dominant factor is the surface 
reference magnetic field B re j . Applying the above con- 
strains, Equation ([7|) reduces to 

B = B ref e-^ z - T] ^2cos(/3x)cos(fy) (8) 

The asymmetric distribution of the non-zero local min- 
ima is periodically spaced, and the magnetic field min- 
ima, B m i n , located at points defined by the coordinates 
{x m in,ymin, dmin), where d min determine the effective z- 
distances of x — y plane that contains the distributed 
lattice sites. Mainly the parameters B m i n and d m i n , de- 
termine the life time of the trapped cold atoms which 
make it critical to locate the d m i n and choose the suit- 
able value of B m i n when micro-fabricating the quantum 
device. The minima locations, along the x, y and z-axis, 
respectively, can be written as 

Xmin = n x a, n x = 0, ±1, ±2, • • ■ (9) 

y m in = n y a(n x ), n y = 0, ±2, ±4, • • • (10) 

dmin « ~ln( B re f + —— ) (11) 

7T y £>z-bias J 

As a result, each individual harmonic potential well 
is localized at its coordinate {x m in,Vmin, dmin) and con- 
fined by magnetic barriers ABt rap where the barriers' 
heights can be estimated to be 

ABj ra = |B^ a J - {BlJ, i = x,y,z (12) 
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TABLE II: Characteristic parameters describe the suitability of the magnetic lattice to host the cold atoms. The magnetic 
walls are produced by the condition T wa u ^ Tbtm- Edge site is denoted by B^ m , while B® jn denotes center site. 



Magnetic lattice parameters 


B x — bias 


R A 


f>B 


d A 


d B 


Simulation output 


G 


0.37 G 


0.5 G 


0.8451 fim 


0.718 ixm 




2 G 


2.03 G 


2.1 G 


0.8451 fim 


0.717 fim 



WW) 

FIG. II. 3: (a) The effect of changing the size of the holes, a^, on the location of the magnetic field local minima along the 
z-axis, B^,i„ located at d m i n above the hole opening center at the surface of the thin film, and (b) the effect of changing the 
periodicity length, a s , across the x — y plane, (c) Simulation result of variating the tunneling barrier heights AB X through 
change in the z-axis magnetization, M z by applying external B z -bias magnetic field in the negative z-direction, and (d) B z -bias 
effect on the gradient of the magnetic sites near the local minima along the z-axis . The thin film thickness: Tbtm = 2 fj,m is 
used in all cases. 



The magnetic fields around the local minima have gra- 
dients in which they are symmetrically distributed across 
the x — y plane. These curvatures, along the x and y axes, 
are given by 
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Due to the x — y symmetry, we find that at the centers 
of the traps, the following holds true 



d 2 B d 2 B 



dx 2 



dy 2 



(14) 



During the trapping process, before atoms reach the 
lowest magnetic field location at the bottom of traps, 
they undergo transverse oscillations. The transverse os- 
cillation frequency v x ^ y depends on the Zeeman sub-levels 
and the curvature of the magnetic field, and it can be ex- 
pressed as 



P / d 2 B k 



k = x,y (15) 



Where gp is the Lande g-factor, fis is Bohr magne- 
ton, F is the atomic hyperfine state with the magnetic 
quantum number mp. The non-zero local minima val- 
ues determine the depth of the harmonic potential wells, 
where the depth A^epth of an individual potential well 
can be expressed as folows 
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Where x = (x,y,z), Kb is the Boltzmann constant, 
^trap an( i Btrap are the maxima and minima values of 
the magnetic field barriers. Increasing the periodicity 
length, a s , between the holes increases the value of the 
non-zero local minima by several Gauss above the Majo- 
rana spin-flip critical minima. Our model shows suitable 
values for a s that positively influence the cold atoms life 
time, where large values of a s causes the magnetic min- 
ima to occur away from the surface, thereby reducing 
the interaction with the surface. The effect of changing 
the periodicity a s , can be used to locate suitable effec- 
tive distances which keep the cold atoms away from the 
Casimir-Polder interacting limits when they are trapped 
near the surface of two-dimensional magnetic lattices, as 
shown in Figure llL37 a-b) . 

The gradient of the trapping magnetic field at each in- 
dividual potential well is of particular importance when 
loading the ultracold atoms into the magnetic lattice. 
Lattice sites with a steeper gradient of their magnetic po- 
tentials maintain a suitable environment with increased 
chances for the destructive thermal Majorana spin-flip 
to exist. Consequently, the loss of ultracold atoms is 
more likely to occur due to the spin-flip process which 
causes an increase in their temperature. In this situa- 
tion the in situ radio frequency (RF) evaporative cooling 
is effectively reducing the number of the re-thermalized 
(hot) atoms and hence enhance the phase space density 
of the compressed MOT to reach the point of BEC tran- 
sition. On the other hand, steep gradients increase the 
chances of having single-atom per lattice site. Moreover, 
the holes opening centers may reduce the effect of the sur- 
face Casimir-Polder interaction when using larger values 
of ah- Simulation results of various dimensions of ah 
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FIG. III.l: (a-d) Size of the magnetic potential wells along the x-axis for different sizes of holes. The simulation is carried out 
using the initial conditions M z = 2 kG, Tbtm = 2 /im and with no external bias fields applied for the different values aj, = a s 
— 1 fim, 3 fim, 5 fim and 7 fim. 

FIG. IV. f: (a) Schematic diagrams of the conditions positive wall, T p ~ wa ii > Tbtm- (b) Density plot representation of the 
simulation results of applying the condition T p ^ wa n — 0.5 fim (j v - wa ii > Tbtm) on the 9x9 magnetic lattice which causes 
reductions in the tilting magnetic potential between the lattice sites, (c) Comparing the magnetic minima locations between 
the center site and the edge site along the z-axis. (d) A fOO x 100 magnetic lattice simulation using the T p - wa u = f fim 
constrain. Random distribution of the sites appears at the edges of the magnetic lattice, (e) Schematic representation of the 
condition negative walls, T n - wa ii < Tbtm, where (f) the density plot representation shows the effect of the condition T n - ma u — 
-0.5 fim on the locations of the sites along the z-axis which are spaced by the tilting potential, (h) For large magnetic lattice 
the effect disappears from the center sites while robustly appears at the edges. As simulation parameters we used M z — 2.8 
kG, Tbtm = 3 fim and no external bias fields has been applied, aj, = a s = 7 fim and = a s — 10 fim is used for large lattice 
size, i.e., n — 100. 



and a s , for a magnetic lattice prepared at initial setup 
are shown in Figure 

Our simulation results show that the unperturbed film 
area, of thickness T wa u, surrounding the n x n matrix 
affect the bottom of most of the potential minima at 
the center of the magnetic lattice. The effect of r wa u 
(magnetic walls effect) on the magnetic lattice sites, in- 
troduces pronounce changes in the value of B min and 
dmin at the edge and center lattice sites. To describe 
the effect of the magnetic walls on some of the character- 
istic lattice parameters, such as the magnetic tunneling 
barrier heights, AB, and the effective distance, d m i n , we 
simulate the trapping magnetic field above the surface of 
the magnetic lattice with different values of T wa u- Con- 
sidering three values for T wa u with respect to the bottom 
layer thickness Tf, tm , the terms positive wall, negative 
wall and surface-equal wall are introduced to denote the 
conditions r p - wa u, T n - Wa ll and T se - wa iu respectively. 

We summarize the two conditions positive wall, 
Tp-waii, and negative wall, T n _ wa n, based on the number 
n of the lattice sites, as follow. Applying the configura- 
tion Tp— W aii results in an unperturbed area with thickness 
that is larger than the thickness of the bottom layer of 
the magn etic la ttice T p _ waU = T waU > n tm , as shown 
in Figure llV.lf a). For small values of n, the T p - wa u 
acts as insitu magnetic bias field sources which surround 
the n x n magnetic lattice. Thus, T v - wa ii configuration 
with proper dimensions may replace the external mag- 
netic bias fields sources created using current carrying 
wires in the BEC experiment. The effect of r p ^ wa u ap- 
pears as reduction in the amount of the tilting magnetic 
potential between sites, as shown in Figure ITV.lf b). Fig- 
ure HEUc) shows the locations of two different magnetic 
minima along z-axis at the center and the edge of the 
magnetic lattice where it is clear the T p - Wa u configura- 
tion reduces the spatial difference between sites. For a 
large magnetic lattice, i.e., n — > 100, the T p - wa ii has no 
effect on the magnetic tunneling barriers and the spatial 



separations between sites at the center of the magnetic 
lattice. However, the edge sites are randomly distributed 
along the z-axis as shown in Figure ITV.lf d). It is also 
noticed that the B m i n of the center sites are close to zero 
and their d m i n s are relatively small but still far from the 
surface. 

When applying the condition T n - W aii = i~ waU < r 6tm , 
the lattice sites are well spaced on the z — x, z — y planes 
as shown in Figure HV.lf e-f). Each m pair of lattice 
sites around the center site have equal B m i n , which led 
us to realize the existence of different magnetic bands 
surrounding the center site. We denote the effect by 
the magnetic band gap structure in a similar way to 
the energy band gap structure in semiconductor devices 
which can be used to simulate interesting problems in 
condensed matter physics using trapped ultracold atoms, 
e.g., ultracold fermions [? ]. Figure HV.lf g'). shows the 
locations of two different sites, i.e., B m i n , along the z- 
axis at the center and the edge of the magnetic lattice. 
The pronounce spatial separation creates the magnetic 
band gap structure in our proposed asymmetrical mag- 
netic lattice. The center of large magnetic lattice is also 
not affected in the case of T n — wa u where center sites have 
BminS close to zero while the edge sites maintain robust 
tilting magnetic potentials. Figure HV.lf h) shows the 
simulation result of applying the condition T n - Wa u on 
a large magnetic lattice with n = 100 sites. 

An asymmetrical magnetic lattice set initially at the 
reference magnetic field, B re f , with no application of ex- 
ternal bias fields, regardless of existence of the wall con- 
strain, creates magnetic confinement with B m i n close to 
zero. This situation can be avoided by changing the peri- 
odicity between the lattice sites, i.e., changing the value 
of the separating period a s between holes to be larger 
than the hole size ah- In Table (|IV[) we summarize the 
simulation results of using different values for the hole 
size, i.e., ah = 0.5 fxm, 1 fim and different values for the 
separating period, a s — 0.5 fim, 1 fim, 1.5 fim, and 2 
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TABLE III: Simulation results show the effect of changing the periodicity a s of a 11 x 11 asymmetric two-dimensional magnetic 
lattice prepared initially at the reference magnetic field B re f and no external magnetic bias fields are applied. 



M z (kG) 


a h (pm) 


a s (um) 










Twaii condition 


3.8 


0.5 


0.5 


0.3 


0.2516 


0.91 


0.2909 
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FIG. IV. 2: The effect of applying the B x -uaa field on the magnetic tunneling barrier AB X , AB Z and the non-zero magnetic 
local minima B^ in (a) along a;-axis and (b) along z-axis for different values of B x -tias- (c) Results of applying external bias 
field B+z-Mas along z-axis, and (d) the effect of applying B~ z ~bias on the magnetic tunneling barriers between the lattice sites, 
where the magnetic lattice is biased by B x ^tiased = 10 G. Simulation input: a s = aj, = 10 am M z = 2.80 kG, Tum = 2 \im 
and T p - wa ii = 1 um. 



Another method of letting B min ^ G, is to apply 
external magnetic bias fields which is commonly used in 
the area of magnetic micro-trapping of cold atoms. Fig- 
ure HV2] shows the effect of applying different magnetic 
bias fields where we simulate the case of ah — a s = 10 
fun. For the tunneling process of the cold atoms to be 
effective periodicity with a relatively small value is de- 
sirable. A detailed analysis of achieving an optimal two- 
dimensional magnetic lattice, using this new method, will 
be reported elsewhere. 

We fabricate the quantum device by milling the n x n 
hole structure using Focused Ion Beam (FIB) on the sur- 
face of magnetic thin film. Figure IIV.3r a-b) show 2D 
image and 3D image of the structure created by using 
Scanning Electron Microscope (SEM) and Atomic Force 
Microscope (AFM), respectively. The the external mag- 
netic field of the asymmetrical magnetic lattice is mea- 
sured using Magnetic Force Microscope (MFM), where 
the structure is biased insitu along the x-axis. Figure 
IIV.3f c-d) show the 3D image of the measured field and 
3D plot of the simulation result, respectively. The exper- 
imental results will be reported elsewhere. 

V. CONCLUSION 

We have proposed a new method to create an asym- 
metric two-dimensional magnetic lattice suitable for 
trapping and confining ultracold atoms and quantum de- 
generate gases prepared in the low magnetic field seeking- 
state. The proposed remanent magnetic structure can be 
fabricated by milling an m x m array of blocks, where 
each block is an array ofnxn square holes. Simulation 
results have shown that it is possible to optimize the two- 
dimensional magnetic lattice, i.e., to achieve symmetri- 
cal distribution of the lattice sites, via the application of 



external magnetic bias fields as well as by applying the 
suitable magnetic walls constrain. The many block fea- 
ture will allow individual addressing for different quan- 
tum states which led us to coin the notion of multi-core 
quantum processors as will be reported elsewhere. We 
have also realized the possibility of creating discrete mag- 
netic levels using the asymmetrical characteristic in the 
distributed magnetic field in which we have denoted the 
discretization effect by the magnetic band gap structure 
on atom chip similar to the energy band gap structure 
in semiconductor devices. This will allow, using the ul- 
tracold quantum degenerate gases, simulating condensed 
matter environment in which many interesting problems 
can be investigated such as exciton formation, Josephson 
effect and long-range entanglement. 

The spacing between the lattice sites is suitable for sin- 
gle atom trapping and detection where well spaced single 
atoms can be trapped at each site. It is of a particular 
importance to highlight the simplicity of the method in- 
troduced in this article for achieving a two-dimensional 
magnetic lattice as it also suits the current state of the art 
technologies. We have shown that applying the external 
bias fields in the negative z-direction controls the tun- 
neling process of cold atoms between sites as well as the 
in-site transitions. Thus, quantum information process- 
ing, such as multipartite entanglement, is also possible 
using such type of magnetic lattices. 
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FIG. IV. 3: (a) Scanning Electron Microscope (SEM) and (b) Atomic Force Microscope (AFM) images of the fabricated two- 
dimensional magnetic lattice, (c) The external magnetic field measured using the Magnetic Force Microscope (MFM) in which 
it agrees with the simulation result in (d) of the insitu biased two-dimensional magnetic lattice. Simulation and experimental 
inputs: a s = aj, = aO /im M z w 2.80 kG, Tbtm = 1 |tm and T p - wa ii = 1 fim. 



